Method and system for cooperative stochastic positioning in a mobile environment

ABSTRACT

Cooperative stochastic positioning in a mobile environment is provided. Position announcements transmitted wirelessly from objects in the immediate area are received at an object of interest or primary object. The position announcements provide the current position data of the respective object in relation to a common coordinate system. The received position announcements are discretizing to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data. Clustering of the data groupings is performed to determine which clusters dataset from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the primary object. A stochastic automata model is then applied to selected cluster dataset to evaluate relative cluster weights in order to determine the relative position of the object of interest. The accuracy of a current position of the object of interest can then be updated based upon determined relative position.

TECHNICAL FIELD

The present disclosure relates to positioning systems and in particular to improving accuracy of absolute and relative positioning between objects in a mobile environment.

BACKGROUND

Regrettably, more than 29,000 fatalities, 2.2M injuries, and $100 Billion dollars in financial losses occur annually on United States roads alone. There has been shared consensus among researchers and governments that those figures can be brought down by applying modern safety applications. The evolving technologies promise to make transportation safer than ever by embedding electronic safety features powered by wireless sensing in vehicles and roads. Location based features are at the heart of this evolution. In particular, identifying the exact position of a moving vehicle (object) is a key aspect to developing most of the vehicular safety features and commercial Location-Based-Applications (LBA).

Global Navigation Satellite Systems (GNSS) such as Geographical Positioning Systems (GPS) have been used for positioning objects with reasonable accuracy. GPS provides typically less than 70 cm accuracy and there have been several studies showing that GPS cannot be efficiently used in urban environment due to dilution of precision and the urban canyon phenomenon. Most safety applications require sub centimeter accuracy with high reliability in urban and sub-urban environments.

Undoubtedly, GPS is the most popular GNSS technology used for positioning objects up-to a normal accuracy of a few meters by timing the transmitted signal along a line-of-sight (LoS) between the satellite and the mobile earth object. If no clear LoS is available between the satellite and the mobile object, ranging to that satellite becomes impossible. The popularity of GPS led to increased interest in Location Based Systems (LBS) where applications behave differently based on user position. Serious and high-end LBS systems such as safety and mission-critical applications cannot tolerate limited positioning accuracy, limited signal availability in urban environment, cloudy/bad weather, or the lack of integrity indicators.

Accurate GNSS positioning (error <30 cm) requires the availability of multiple satellite signals (5+) which is impossible in urban and metropolitan areas. The transportation industry has an inevitable and imperative need to resolve the problem of inaccurate positioning in order to unlock the essential development of safety and automation applications. Unfortunately, GNSS and GPS systems exhibit the following inherent distinctive limitations:

a) Limited signal availability in urban environments, cloudy skies, or tunnels.

b) Insufficient accuracy to serve serious and high-end LBS systems like safety and mission-critical application.

c) Loss of precision until Time-To-First-Fix (TTFF) is made available. The TTFF is the time required for the receiver to acquire ephemeris as well as an almanac for all satellites that contains coarse orbit and status information for each satellite in the constellation. Fix times are unacceptably long, and fixes may never be reached when attenuation levels exceed 30 dB, which is likely in an urban environment and in bad weather.

d) Limited redundancy since GPS has no alternative system to be used in its absence.

e) Satellite-based systems are too centralized and lack the desired localized control.

Accordingly, there is a need for an accurate positioning system and method that provides improved positioning accuracy in a mobile environment.

BRIEF DESCRIPTION OF THE DRAWINGS

Further features and advantages of the present disclosure will become apparent from the following detailed description, taken in combination with the appended drawings, in which:

FIG. 1 shows a representation of mobile environment with random motion where (a) shows before, and (b) shows after the motion;

FIG. 2 shows a method for cooperative stochastic positioning in a mobile environment;

FIG. 3 shows a detailed method for cooperative stochastic positioning in a mobile environment;

FIG. 4 shows a representation of an example of multi-path propagation channel model;

FIG. 5 shows method for performing Advanced Fuzzy Clustering Algorithm;

FIG. 6 shows system for performing cooperative stochastic positioning in a vehicular mobile environment utilizing dedicated short-range communications (DSRC); and

FIG. 7 shows system for performing cooperative stochastic positioning in a mobile environment between objects.

It will be noted that throughout the appended drawings, like features are identified by like reference numerals.

DETAILED DESCRIPTION

Embodiments are described below, by way of example only, with reference to FIGS. 1 to 7.

In accordance with an aspect of the present disclosure there is provided a method for cooperative stochastic positioning in a mobile environment executed by a processor in an object of interest, the method comprising: receiving position announcements transmitted wirelessly from a plurality of objects, the position announcements providing the current position data of the respective object in relation to a common coordinate system; discretizing the received position announcements to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data; performing clustering of the data groupings to determine which cluster datasets from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the object of interest; applying stochastic automata model to selected cluster dataset to evaluate relative cluster weights in order to determine the relative position of the object of interest; and updating accuracy of a current position of the object of interest based upon determined relative position.

In accordance with another aspect of the present disclosure there is provided a system for cooperative stochastic positioning in a mobile environment the system comprising: a first receiver for receiving positioning announcements from a plurality of objects, the position announcements providing the current position data of the respective object in relation to a common coordinate system; a processor coupled to the first receiver; and a memory comprising instructions for execution by the processor, the instructions comprising: discretizing the received position announcements to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data; performing clustering of the data groupings to determine which cluster datasets from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the object of interest; applying stochastic automata model to selected cluster dataset to evaluate relative cluster weights in order to determine the relative position of an object of interest; and updating accuracy of a current position of the object of interest based upon determined relative position.

In accordance with yet another aspect of the present disclosure there is provided A computer readable memory providing instructions for performing cooperative stochastic positioning in a mobile environment, the when executed by a processor in an object of interest performing the method comprising: receiving wirelessly, position announcements transmitted from a plurality of objects, the position announcements providing the current position data of the respective object in relation to a common coordinate system; discretizing the received position announcements to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data; performing clustering of the data groupings to determine which cluster datasets from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the object of interest; applying stochastic automata model to selected cluster dataset position to evaluate relative cluster weights in order to determine the relative position of the object of interest; and updating accuracy of a current position of the object of interest based upon determined relative position.

The present disclosure provides a method and system for accurate positioning in a mobile environment using positioning information cooperatively from objects in the relative area. The disclosure is applicable to implementations of advanced safety vehicular applications but also equally applicable to any environment where objects are mobile and position data is required. In response to the immanent need for higher accuracy and integrity, a simple stochastic approach is provided that can improve and enhance vehicular or object positioning.

Suppose in an environment, objects are moving freely in any direction in space. At any point the position or location of all objects, as described in FIG. 1 a objects moving within an environment are shown at an instance in time 110. At a second instance in time 120 the position of the objects can be at different positions as described in FIG. 1 b. At any given point in time each object is located in a unique position defined by three coordinates Λ_(i)=(x_(i), y_(i), z_(i)) relative to a known absolute (or relative) position. It should be noted that the selected set of coordinates could be angular, Cartesian, or any other coordinate system. The coordinate system may be defined based upon the type of object and the environment in which positions determination is required. For example, vehicles may use GPS where as packages may use coordinate system defined relative to a building structure.

Discrete Time Announcements

Assume that many objects are announcing, via a wireless technology, their positions to the best of their knowledge to be Λ={Λ₁, Λ₂, . . . , Λ_(n)}. Then assume that the Object-of-Interest (OOI), as shown in FIG. 1, either knows its position, but is interested in utilizing the available information to improve its accuracy, or completely unaware of its position. The disclosed system and method provides the OOI with the ability to better define, or improve knowledge about, its position. This can be achieved by finding Λ_(OOI) at time (t) and timing accuracy (ζ), and with the highest positioning accuracy possible (minimum ΔΛ_(min)).

This would be quite easy if OOI can calculate the distance (distances between objects are typically called range) between itself and each of the objects announcing their positions. In other words, the set of ranges

={

,

, . . . ,

} define the range of OOI to each of the announcing objects OX₁ to OX₇. Therefore, if Λ and ρ are both available at time (t) (timing accuracy ζ), a simple multi-lateration would be possible and a simple Euclidean math would lead us to the solution. Even if the distances between all objects are too long, corrections to the Euclidean math are known to fix those errors. For clarity ρ_(i) is the projection of the range between OOI and object OXi on the set of chosen coordinates and can be written as

_(i)=(ρx_(i), ρy_(i), ρz_(i)). Now a typical mathematical approach can be simplified as follows:

Knowing the sets:

Λ={Λ₁,Λ₂, . . . ,Λ_(n)} &

={

₁,

₂, . . .

_(n)};  (1)

The exact coordinates of OOI (Λ_(OOI)=(x_(OOI), y_(OOI), z_(OOI))) can be obtained as

x _(OOI) =f(Λ,

) y _(OOI) =f(Λ,

) z _(OOI) =f(Λ,

)  (2)

Unfortunately, a set of inherent practical limitations lead to the following:

-   -   i. The correct definition of Λ_(i) would be:         -   Λ_(i)=(x_(i)+δx_(i), y_(i)+δ_(i), z_(i)+δz_(i)) where             (δx_(i), δy_(i), and δz_(u)) are the projected errors in             estimated position of object i (OOI), and that each δ is a             range (positive to negative). The error (δ) is an indication             of the accuracy of equipment used and the time since last             accurate position was measured. Errors in Λ_(i) may or may             not be projected for each coordinate. Yet, since error             values are relatively small compared to (x_(i), y_(i),             z_(i)) and the range (ρx_(i), ρy_(i), ρz_(i)), if the error             in estimated position is not available for each coordinate,             a simple approximation can be obtained by assuming the             Euclidian vector of the error.     -   ii. Assuming that the OOI can build some sort of confidence in         the announced position Λ_(i) received from a particular object,         lets' call that εΛ_(i). Now it is important to distinguish         εΛ_(i) from δΛ_(i). For one, the value of εΛ_(i) represent the         infinite range (0-1) where zero indicates no confidence and one         indicates complete confidence. By that, εΛ_(i) indicates a level         of trust in the received information. Should an object attempt         to deceive a group of objects to build wrong conclusion on their         position, εΛ_(i) can capture that in different ways. As an         example, OOI may develop a lower confidence value in Λ_(i).         (εΛ_(i)<9.0 e 20) if the object i (OOI) position information         arrived missing security credentials. Further, an object that         continues to announce position information that reflect lower         information quality from other objects will get diminishing         value for εΛ_(i). Alternatively, objects that continue to         announce information that fits the collective sense gathered         compared to other objects will get increasing εΛ_(i).     -   iii. During the calculation of the range ρi, each range should         have been presented as:         -   _(i)=(ρx_(i)+δρx_(i), rρ_(i)+δρy_(i), ρz_(i)+δρz_(i)) where             (δρx_(i), δρy_(i), and δρz_(i)) are the projected errors in             estimated range of object i (OOI), and that each δ is a             range (positive to negative). The error is an indication of             the accuracy of equipment used. Errors in ρi may or may not             be projected for each coordinate. Yet, since error values             are relatively small compared to (x_(i), y_(i), z_(i)) and             the range (ρx_(i), ρy_(i), ρz_(i)), if the error in             estimated range is not available for each coordinate, a             simple approximation can be used by assuming the Euclidian             vector of the error.

The stochastic solution for positioning problem has the following characteristics:

-   -   1) If the number of objects announcing their position (n) is too         low, the accuracy of the stochastic approach would suffer.         Similarly, if (n) is too high, it is better to focus the         stochastic solution on the best trusted, closer objects with         minimum chances of errors. The system and method includes         built-in measures to monitor the upper and lower limits for (n).     -   2) The distribution of the confidence in a particular measure         εΛ_(i) is dynamic, and follows a feed forward approach. This         requires maintaining a history of previous nodes accuracy, and         therefore, makes the algorithm memory demanding.

3) The method design strongly guards against attacks; at least (½ n) must provide (consistent) malicious information before the OOI would show some signs of lower accuracy. Even in that case, it is fairly simple to predict malicious objects and isolate them. In other words, it is extremely difficult to mislead the method by faking invalid position messages.

-   -   4) By just listening to exchanged messages within the         environment, the OOI can draw pretty accurate image about its         existence.

As illustrated, the system and method provides higher accuracy by merely utilizing available information. It may require fairly large memory storage, but has been proven to show strong resilience to malicious attacks.

FIG. 2 shows a method for cooperative stochastic positioning in a mobile environment. The method is implemented in a receiving or primary object to aid in improving the accuracy of the location of the object. The object of interest receives position information (202) wirelessly for other object in the immediate area. The position information defines each objects current location data relative to a standard coordinate system, such as for example GPS. The received data can be discretized into a number of clusters (204) based upon a chosen criteria such as for example time varying frequency relative to the arrival times of the position information. Ranging information can also be simultaneously determined (206) and associated with the received position information. The ranging information can be determined by the receiver of the positioning information or simultaneously determined by other ranging technologies and correlated with the received position information. Alternatively the ranging information can be generated by applying a channel propagation model to the received position information. Data clustering is then performed (208) on each cluster to classify the clusters into competing data sets based upon the chosen criteria. The criteria may also be associated with previous accuracy of location information provided by an object. For example objects that have previously provided low accuracy data may be clustered together or weighted differently in the cluster datasets. Clusters that provide sufficient accuracy within a defined distance can than be processed using a stochastic automata model (210) to perform multi-lateration to determine a more accurate position of the object of interest. The objects position can then be updated based upon the determined position (212) and the process continues at (202) for the next set of received position data.

FIG. 3 shows a detailed method for cooperative stochastic positioning in a mobile environment. The method is implemented in a receiving object to aid in improving the accuracy of the position of the object. The object receives position information (202) for other object in the immediate area. The position information defines each objects current position data relative to a standard coordinate system, such as for example GPS. The received information contains time and position information (302), however the information can be possibly inaccurate information or unavailable altogether (304). The received data can be discretized into a number of groupings (204) based upon time varying frequency relative to the arrival times of the position information. Ranging information can also be simultaneously determined (206) and associated with the received position information. The ranging information can be determined by the receiver of primary object or simultaneously determined by other ranging technologies and correlated with the received position information. Alternatively the ranging information can be generated by applying a channel propagation model to the received position information. Data clustering is then performed (208) on each cluster to classify the clusters into competing data sets. Classifiers are applied to the determined clusters (306) based on the selected clustering criteria. For each cluster, if the accuracy of the cluster is less than the required accuracy outside a defined threshold, (NO at 308), the cluster can be divided (310) and classifiers applied (306). For clusters that provide sufficient accuracy within a defined threshold, (YES at 308), they can than be processed using a stochastic automata model (210) to perform multi-lateration to determine a more accurate position. A stochastic automata model can then be applied (312) to the clusters to identify stronger data weights in the position data. The weighted multi-lateration position can be calculated using the weighted position data (314). The calculated position data can then be validated (316) against the existing position data of the primary object. If the calculated position data is less than an accuracy threshold (<required accuracy threshold at 318), the least weighted data set can be removed from the data (318). If the data provides sufficient accuracy (>=required accuracy at 318) the objects position can then be updated (212) and the process continues at (202) for the next set of received position data.

Continuous Time Announcements

It would be ideal if the OOI receives position announcements from all (n) objects at discrete time points. Unfortunately, this is impossible. Objects are expected to announce their positions sporadically and randomly. Therefore, the method discretizes the available information at a time varying frequency. Since discretizing available information would involve extrapolation, it can affect accuracy. The effect on position accuracy can be accommodated by updating δ values. The method must also make a decision on the best time to recalculate position.

Assuming the Last-Known-Accurate-Position (LKAP) is Λ₀=(x₀, y₀, z₀), and given the position expectation at the next point Λ₁=(x₁, y₁, z₁), the OOI would discretize available information if the time since Λ₀ exceeds mζ, or if the estimated errors in Λ₁ contradicting available information exceeds acceptable tolerance. This system and method provides ways to mitigate the timing granularity effect mc and to decide when to re-evaluate the new position Λ.

Channel Propagation Model

Assume a wireless signal has N resolvable propagation paths between the transmitter and receiver. Since each reflector propagates multiple signals, and there are potentially have multiple reflected signals, the strongest reflected signal is selected. The multipath signal parameters used are: Angle of Departure (AoD φ_(i)), Angle of Arrival (AoA θ_(i)), and delay of Arrival (DoA τ_(i)). All of (φ_(i), θ_(i), and τ_(i)) can be measured using any available technique with respect to common bearing direction. Then, as illustrated in FIG. 4, let (x_(Λ), y_(Λ)), (x_(o), y_(o)), and (x_(i), y_(i)) be the true position of, respectively, the signal source 410, the object of interest 420, and the ith reflective point where (x_(o), y_(o)), and (x_(i), y_(i)) are unknown. Let ρ′_(i), p″_(i) be the lengths of the segments forming the ith path respectively. Finally, let φ_(i) and θ_(i) be the angles of departure and arrival for the ith path from the source 410 to the object of interest 420.

From FIG. 4, φ_(i) and θ_(i) are obtained as a function of the locations of the mobile target and the reflection point.

$\begin{matrix} {{{\theta_{i}\left( {x_{o},y_{o},x_{i},y_{i}} \right)} = {\arctan \left( \frac{y_{i} - y_{o}}{x_{i} - x_{o}} \right)}}{{\phi_{i}\left( {x_{o},y_{o},x_{i},y_{i}} \right)} = {\arctan \left( \frac{y_{i} - y_{s}}{x_{i} - x_{s}} \right)}}} & (3) \end{matrix}$

Equation 3 applies for 1=1, . . . , N. Then assuming c is the, corrected, propagation speed, the Time Difference of Arrival (TDoA) can be determined by:

$\begin{matrix} {{{\tau_{i}\left( {x_{o},y_{o},x_{i},y_{i}} \right)} = \left( \frac{\rho_{i} - \rho_{1}}{c} \right)}{and}{{i = 2},\ldots \mspace{14mu},N}} & (4) \end{matrix}$

where ρ_(i)=ρ′_(i)+ρ″_(i), and:

ρ′_(i)=√{square root over ((x _(i) −x _(o))₂+(y _(i) −y _(o))²)}{square root over ((x _(i) −x _(o))₂+(y _(i) −y _(o))²)}

ρ″_(i)=√{square root over ((x _(i) −x _(s))²+(y _(i) −y _(s))²)}{square root over ((x _(i) −x _(s))²+(y _(i) −y _(s))²)}  (5)

Since the unknown position (x_(o), y_(o)) needs to be obtained from the known position (x_(Λ) y_(Λ)), given the uncertainty in ({circumflex over (φ)}_(i), {circumflex over (θ)}_(i), and {circumflex over (τ)}_(i)), The expected statistical error in measuring (φ_(i), θ_(i) and τ_(i)) is applied such that:

{circumflex over (τ)}_(i)=τ_(i)(x _(o) ,y _(o) ,x _(i) ,y _(i))+n _(τ) _(i)

{circumflex over (θ)}_(i)=θ_(i)(x _(o) ,y _(o) ,x _(i) ,y _(i))+n _(θ) _(i)

{circumflex over (φ)}_(i)=φ_(i)(x _(o) ,y _(o) ,x _(i) ,y _(i))+n _(φ) _(i)   (6)

where (n_(φ) _(i) , n_(θ) _(i) , and n_(τ) _(i) ) are the statistical errors, i=1, . . . , N for {circumflex over (φ)}_(i) and {circumflex over (θ)}_(i), but i=2, . . . , N {circumflex over (τ)}_(i).

Therefore, when the number of paths N≧3, (3N−1) measurements and (2N+2) unknown parameters are determined. The problem yields a non-linear estimation problem that can be solved using machine learning or stochastic learning automata. Using a stochastic learning automata all nodes cooperate to arrive to better relative ranging. Therefore, the collective behavior of the selected set of nodes, or a cluster, are guaranteed to converge, since the accumulation of relative distances will naturally tend to marginalize low accuracy ranging in favor of the multiplicity of better accuracy ranging. Further, since the boundaries of the selected set of nodes (a cluster of nodes) are finite, the problem lends itself to deterministic finite automata. In the following section the formulation of the automata model is disclosed.

The approach used in this section defines a way to handle wireless channel propagation model that is efficient for many wireless technologies, but remain isolated from the stochastic method. Manufacturers might use alternate approaches to calculate the range between OOI and OX. The method that is illustrated here is given to show the use of non-Cartesian coordinates with the stochastic approach. Alternative approaches can be used for instance, using an infrared ranging method combined with the disclosed stochastic approach that utilizes the calculated ranges but not to the ranging method in itself.

Clustering

A soft classifier is utilized to split the set of input values A into competing sets or clusters. Each set of Λ would include closer consistency in its parameters. In other words, elements of a set Λ_(a) when combined together, they end up with a crisp value for Λ_(OOI) that has the following features:

-   -   1—All elements of Λ_(a) can be used to find Λ_(OOI) with         accuracy <=ΔΛ at time (t) and timing accuracy (ζ).     -   2—Elements that violate the previous condition may join Λ_(b),         where all elements of Λ_(b) can be used to find Λ_(OOI) with         accuracy <=ΔΛ at time (t) and timing accuracy (ζ). Same applies         to compute Λ_(c), Λ_(d), . . . and so on.     -   3—At the end, the remaining a set of elements will represent         elements that would together violate the required accuracy         measures. Those elements are ignored or used for swapping only.         The above mentioned function can be achieved by applying fuzzy         classifiers, rough sets, neural networks, or other similar         classifiers. As disclosed in FIG. 5, one classifier method is an         Advanced Fuzzy Clustering Algorithm (AFCA) that performs the         required classification to performing data clustering (308).

The method commences with the selection of clustering criteria (502) such as least-mean-square (LMS), arrival delay, data confidence or any other grouping criteria that can be utilized inferable from the received position data. The clustering criteria may also be associated with previously data received from an object to determine an accuracy of the data. In such a case a confidence interval may be assigned to position information or other characteristics to define accuracy of the data. Based upon the configuration of the system, the criteria may be fixed, for example during startup of the system or may provide multiple criteria options that are selected when a particular criteria for clustering does not provide accurate results and the method reinitiated, The cluster data set are selected (504) from the data grouping based upon the chosen criteria. The selected criteria is then used to minimize error in cluster accuracy (506) by switching one data position within the cluster at a time. This can be performed for example by, letting the given data set consist of N data points, denoted by the vectors:

$\begin{matrix} {X_{i} = \begin{bmatrix} x_{i\; 1} \\ x_{i\; 2} \\ \vdots \\ \vdots \\ x_{ip} \end{bmatrix}} & (7) \end{matrix}$

Where i=1, 2, . . . , N. The AFCA classifies the data into R disjoint groups θ₁, θ₂, . . . , θ_(R), according to a certain selected criterion. In order to elaborate this point, assuming the center of gravity of the total vectors being at the origin, i.e., equation 3 can be converted to:

$\begin{matrix} {\overset{\_}{X} = {{\frac{1}{N}{\sum\limits_{i = 1}^{N}X_{i}}} = 0}} & (8) \end{matrix}$

Otherwise, a transformation from X_(i) to X requires that the transformed pattern satisfies equation 8, and the normalized total scatter matrix T based on θ is given by:

$\begin{matrix} {T = {\frac{1}{N}{\sum\limits_{i = 1}^{N}{X_{i}X_{i}^{t}}}}} & (9) \end{matrix}$

Where the superscript t, indicates a transposed matrix, and the normalized intra-group scatter for class θ_(j) is:

$\begin{matrix} {\Gamma_{j} = {\frac{1}{N_{j}}{\sum\limits_{x_{i} \in \theta_{j}}\; {\left( {X_{i} - {\overset{\_}{X}}_{j}} \right)\left( {X_{i} - {\overset{\_}{X}}_{j}} \right)^{t}}}}} & (10) \end{matrix}$

Where N_(j) is the number of vectors in group θ_(j), and:

$\begin{matrix} {{\overset{\_}{X}}_{j} = {\frac{1}{N_{j}}{\sum\limits_{x_{i} \in \theta_{j}}\; X_{i}}}} & (11) \end{matrix}$

Therefore, the average normalized intra-group scatter (508) for the total data set θ is:

$\begin{matrix} {\Gamma = {\sum\limits_{j \approx 1}^{R}\; {p_{j}\Gamma_{j}}}} & (12) \end{matrix}$

Where p_(j)=N_(j)/N and the normalized intra-group scatter (510) is given by:

$\begin{matrix} {B = {\sum\limits_{j = 1}^{R}\; {p_{j}{\overset{\_}{X}}_{j}{\overset{\_}{X}}_{j}^{t}}}} & (13) \end{matrix}$

Therefore, the following matrix identity is defined:

T=Γ+B  (14)

Since different data criterion that can be expressed in terms of matrices T, Γ, and B, a chosen criteria to minimize is:

$\begin{matrix} {S_{0} = {{{tr}\; \Gamma} = {\sum\limits_{j = 1}^{R}\; {\sum\limits_{x_{i} \in \theta_{j}}\; {{X_{i} - {\overset{\_}{X}}_{j}}}^{2}}}}} & (15) \end{matrix}$

Where ∥ . . . ∥ is the Euclidean norm and tr Γ is the trace of Γ. It is noted that S₀ is only invariant under an orthogonal transformation of the data space, and the AFCA method based on S₀ is most appropriate for fairly concentrated clusters. The two clustering criteria which are invariant under any linear non-singular transformation are to maximize:

$\begin{matrix} {S_{1} = {\frac{\det \; T}{\det \; \Gamma} = {\sum\limits_{k = 1}^{p}\; \left( {1 + \lambda_{k}} \right)}}} & (16) \\ {S_{2} = {{\det \; \Gamma^{- 1}B} = {\sum\limits_{k \approx 1}^{p}\; \lambda_{k}}}} & (17) \end{matrix}$

Where λ_(k), k=1, 2, . . . p are the eigenvalues of the matrix Γ⁻¹ B. However, the AFCA clustering algorithm is based on the three criteria and cannot have reasonable performance and computational complexity simultaneously. It is easy to show that under a normalizing transformation of the data space, any alternative criteria is suboptimal, i.e. when

X→Z=AX  (18)

Where A is a linear non-singular transformation matrix such that

ATA ^(t) =I  (19)

The elaborated criterion demonstrate that the criterion S′₀, which is the criterion S₀ after the normalizing transformation of equation 18 is applied, has the advantages of both optimum performance and low computational complexity compared to the criteria S₀, S₁, and S₂.

For presentational purposes only, the least-mean-square (LMS) criterion is used. The main advantages of this criterion are that the method is invariant under any linear non-singular transformation and classification can be done by either a linear or a generalized linear machine. Yet it is important to indicate that the AFCA may use alternative S′₀ criterion especially when clusters are not close to equi-probable like the case with the short-range signal. Now let:

Φ(X)=[Φ₁(X),Φ₂ (X), . . . ,Φ_(q)(X)]^(t)  (20)

Where Φ_(i)(X), i=1, 2, . . . , q, are linearly independent, real, single-valued and continuous functions of the components of X. The corresponding Euclidean space of Φ(X), is E^(q). The LMS criterion is to minimize the quantity:

$\begin{matrix} {{h(M)} = {\frac{1}{2}{\sum\limits_{j \approx 1}^{R}\; {p_{j}E_{j}{{{M\; {\Phi (X)}} - \alpha_{j}}}^{2}}}}} & (21) \end{matrix}$

Where M is an (R−1)×q weight matrix, α_(j), j=1, 2, . . . R, are the reference points with the properties that:

$\begin{matrix} {{\alpha_{j}} = 1} & (22) \\ {{{\alpha_{j}^{t}\alpha_{k}} = \frac{1}{R - 1}},{j \neq k}} & (23) \\ {{{{{\alpha_{j} - \alpha_{k}}} = {{\alpha_{j} - \alpha_{i}}}},{and}}{{{i \neq j \neq k} = 1},2,\ldots \mspace{14mu},R}} & (24) \end{matrix}$

And E_(j) is the expectation taken with the probability distribution of X's in group θ_(j). The selected criteria can then be assessed (512) to determine if it should be tightened/loosened or if additional clustering criteria should be selected with the method restarting at (502) to select new clustering criteria, (204) in FIG. 2, if a new criteria is to be used. It can then be determined if clusters can be combined (514) if the data provided by each cluster is within a defined position accuracy relative to the defined coordinate system, for example 20 mm. If clusters can be combined, (YES at 514), the clusters are both within the desired accuracy, the selected criteria can be used to minimize error of the combined cluster (506). If the clusters can not be combined, (NO at 514), with each other based upon relative position accuracy, it is then determined if the cluster accuracy is within the required accuracy (516) of the system. If the cluster is within the required accuracy, for example 20 mm of the known position, (YES at 516), the cluster can then be provided to the stochastic automata (518). If the accuracy is not within the required accuracy (NO at 516), it is then determined if the cluster can be divided (520) based upon the relative accuracy of the cluster data. If the cluster can be divided (YES at 520), the divided cluster is provided to select the dataset and clustering criteria (504). If the cluster cannot be divided (NO at 520), the cluster is removed (524). This process is repeated until all the data is analyzed, the position data either being included in a cluster that is passed to the stochastic automata or removed from the dataset.

The following sections describe the core stochastic method for calculating Λ.

Formulation of Automata Model

A classical deterministic finite automata problem is summarized, and then, the stochastic automaton is formulated that is considered in mobile environment situations. Alternative approaches can be used to calculate the range or classify calculated ranging information. Yet, the stochastic approach defined in this section and the following sections apply on how to calculate the high accuracy position of OOI. A zero sum is desired where all cooperative nodes (or participants) continuously reinforce lower square law. The next round of readings to restart the competition. Therefore, there will be a way to stop the competition or to indicate the point of departure.

The finite deterministic automaton is defined by a quintuple {Λ,

, φ, g, h}, where Λ is the set of position inputs,

is the set of estimated ranges, and φ is the set of confidence values where the higher the confidence the better. The set φ is treated as a discrete set of states by enforcing disconnected steps of 0.10, and therefore, maintains a finite state. The corresponding lowercase letters denotes members of the defined sets (e.g., λ_(i), ρ_(i), and φ_(α)). Now g is the output function, and therefore, ρ(t)=g[φ(t)]. Similarly, the next-state function φ(t+1)=h[φ(t), λ(t+1)]. The next equation is the canonical equation for finite automaton. Restricting the input set to two values, 0 and 1 respectively, called non-penalty and penalty where the penalty applies if the position exceeds certain threshold from the expected position. Therefore, I set λ₁=0 and λ₂=1. Finally, the random medium can be described by: C=C(p₁, p₂, . . . , p_(k)), where:

p _(α) =Pr[λ(t+1)=λ₂|ρ(t)=ρ_(α)]  (25)

Since the input to the automaton at instant (t+1) is the output of the medium at instant t. When p_(α)s are constant, a stationary random state is defined. The state function h determines the learning behavior of the automaton. With suitable choice of h, the average penalty of the automaton decreases with time to an asymptotic value.

The stochastic automaton is defined by the sextuple {Λ,

, φ, g, Π, T}. Similarly, Λ would be the set of two inputs (1 and 0:1 for penalty and 0 for non-penalty),

is the set of estimated ρ ranges, φ is the set of states, g is the output function ρ(t)=g[φ(t)] which is one-to-one deterministic mapping from state set to the output set. Now, since the number of outputs equals the number of states, the vector Π is the state probability vector and Π(t)=(Π₁(t), . . . ,Π₁(t)) controlling the choice of the state and hence the output at instant t. Therefore, the state φ_(i) would be chosen at instant t with probability:

$\begin{matrix} {{\prod\limits_{i}\; (t)},{{\sum\limits_{i = 1}^{r}\; {\prod\limits_{i}\; (t)}} = 1}} & (26) \end{matrix}$

Finally, T defines the reinforcement scheme which drives Π(t+1) from Π(t). This can be simplified as:

Π(t+1)=T[λ(t),φ(t)]Π(t)  (27)

Where T is not an explicit function of t, but t is the discrete-time parameter. Depending on the previous action and the environmental response to the action taken, the reinforcement scheme changes the probabilities controlling the choice of the state at the next-round. With suitably designed reinforcement scheme, the average penalty tends to decrease with time. This guarantees a conversion with relative ranging in a stationary state. However, in a dynamic environment, the timing may not be enough to reach sufficiently low average penalty. This can be mitigated by decreasing the time granularity to increase the number of rounds. Computational cost here is pretty minimal, and it is considered a good extension for this research.

In order to perform the linear and nonlinear reinforcement you need to consider the effect of environmental mobility on the performance of the proposed stochastic automaton. This resembles pretty much the operations of adaptive controller in a mobile environment. The stochastic state probabilities are continuously altered according to a reinforcement scheme in response to penalties received from the environment. The automaton adapts by reducing the average penalty. Periodic perturbations of penalty strengths are used as test-singles to drive analytic expressions describing the tracking behavior or the automaton operating under a linear reinforcement scheme. It can be proven that the parameters that determine the ability of a stochastic automaton to track the, perceived, mobility of the environment are the Eigen-values of the transition matrix R(α, β) as illustrated later.

Therefore, changes in the environment are reflected as perturbations in the asymptotic state probability values of the automaton. In this manner, the automaton is said to track the environment. However, as the perturbation frequency (ε) increases to such an extent that the automaton starts responding to the average value of perturbed parameter, then the automaton loses the ability to track. Hence, it is important to express the upper limit of the perturbation frequency (ε_(u)) which keeps the automaton tracking the environment mobility. These limits are expressed as functions of the Eigen-values of the transition matrix R(α, β). The Eigen-values also determine the correlation between any two state probability values separated by a number of transition intervals. The smaller the correlation, the higher is the upper limit (ε_(u)). For the purpose only on the square law reinforcement scheme is discussed which can be defined as follows. If φ(t)=φ_(i), in other words if ρ(t)=ρ_(i), then:

$\begin{matrix} \begin{matrix} {{{{\prod\limits_{i}\; \left( {t + 1} \right)} = {\prod\limits_{i}^{2}\; (t)}},}\;} & {{{if}\mspace{14mu} \lambda} = 1} \\ {{{\prod\limits_{i}\; \left( {t + 1} \right)} = {1 - {\sum\limits_{j \neq i}\; {\prod\limits_{j}^{2}\; (t)}}}},} & {{{if}\mspace{14mu} \lambda} = 0} \end{matrix} & (28) \\ {{{{and}\mspace{14mu} {\forall{j \neq i}}};}\begin{matrix} {{{\prod\limits_{j}\; \left( {t + 1} \right)} = {{\prod\limits_{j}\; (t)} + \frac{{\prod\limits_{i}\; (t)} - {\prod\limits_{i}^{2}\; (t)}}{r - 1}}},} & {{{if}\mspace{14mu} \lambda} = 1} \\ {{{\prod\limits_{j}\; \left( {t + 1} \right)} = {\prod\limits_{j}^{2}\; (t)}},} & {{{if}\mspace{14mu} \lambda} = 0} \end{matrix}} & (29) \end{matrix}$

Using this reinforcement, the steady-state performance of the automaton can be summarized as follows:

$\begin{matrix} {{{{if}\mspace{14mu} p_{i}} < \frac{1}{2}},{p_{j} > \frac{1}{2}},{\therefore{j \neq i}}} & (30) \end{matrix}$

Also the vector Π^((i)), which is defined by Π_(i) ^((i))=1 and Π_(j) ^((i))=0, where j≠i, would control the choice of states at the steady state. Therefore, the automaton actually chooses the probability of the state corresponding to the lowest penalty probability. This is promoted as the best possible performance. In this case, the automaton will be said to perform optimally, given that the state is assumed stationary at each value of t. If however, all p_(i)>½, the stable state probability vector defining the choice of states asymptotically. Then, by reapplying Eigen-values on the transition matrix as performed before, and following the same approach to a great degree of accuracy by Π*, where:

$\begin{matrix} {\prod\limits_{i}^{*}\; {= \frac{\frac{1}{{2\; p_{i}} - 1}}{\sum\; \frac{1}{{2\; p_{j}} - 1}}}} & (31) \end{matrix}$

Equation 31 is valid when the number of states is two. From the same equation it can be derived that the state corresponding to a lower penalty probability has higher probability of being chosen. This performance is referred to as expedient. In the case when more than one p_(i)<½, it leads to a situation where suboptimal performance is possible and might converge, it can only happen if a major malicious attack has been attempted successfully on too many objects. However, it is highly unlikely that so many objects can be led to invalid p_(i) unless the entire system stops for a long period that is much higher than (ε_(u)). Since this situation is practically unachievable, Π* will always converge.

Exploration of a Scenario

Now to run through a case showing how the competition process works, some definitions are introduced. For instance, the penalty of the preceding round can be called unit-loss while the non-penalty would be called unit-win. The ith output is identified and associated with the ith estimated range, therefore, p_(i) is the probability of unit-loss based on range ρ_(i), and q_(i)=1−p_(i) is the probability of unit win for range ρ_(i). In that sense, the winning node would define the correct absolute-coordinates. Further, winning ranges would define the second best set of relatively well-positioned nodes. Therefore, relatively well-positioned nodes can adjust their position estimate based on the winning node absolute-coordinates and winning ranges by using multi-lateration. Finally, badly positioned nodes are those nodes with subsequent unit-loss. They can fix their coordinates based on winning ranges from either the winning node or the well positioned nodes.

Let us consider two automata, Σ¹ and Σ² taking part in the competition process. The input λ^(j)(t) of the automaton Σ^(j) has two values λ₁ ^(j)=0, i.e. a win, and λ₂ ^(j)=1, i.e., a loss. If ρ_(j) is the number of states for Σ^(j), then, for the stochastic automata defined in the previous section, the number of ranges is also ρ_(j), (i.e. ranges are ρ₁ ^((j)), ρ₂ ^((j)), . . . , ρ_(rj) ^((j)). Hence, if

(t) at time t is the set {ρ⁽¹⁾(t), ρ⁽²⁾(t)}. The outcome of the round at time t would be the set Λ(t+1)={λ⁽¹⁾(t+1), λ⁽²⁾(t+1)}. A competition Γ with automata Σ¹ and Σ² is defined if, for all sets

(t), the probabilities P[

(t), Λ(t+1)] of its outcomes Λ(t+1) are given. In the usual meaning of the competition theory, the payoff functions M^(j)(

), j=1, 2, defining a competition Γ*, denote the expectation of winnings for the jth node, the set of ranges being denoted by F, and are related to P(

, Λ) in the following way:

$\begin{matrix} {{{M^{j}(F)} = {\sum\limits_{\lambda^{(i)}}\; \left\{ {{P\left\lbrack {{F;\lambda^{(i)}},0} \right\rbrack} - {P\left\lbrack {F,\lambda^{(i)},1} \right\rbrack}} \right\}}},{i \neq j},} & (32) \end{matrix}$

In a zero sum situation:

$\begin{matrix} {{\sum\limits_{j}\; {M^{j}(F)}} = 0} & (33) \end{matrix}$

If at time t both automata Σ¹ and Σ² are in states φ_(α) ⁽¹⁾ and φ_(β) ⁽²⁾, the system is in state (α, β). If the R(α, β) represents the final probability of the system state being (α, β), then W^(i), the expected value of the winning of Σ^(j) is given by:

$\begin{matrix} {{W^{i} = {\sum\limits_{\alpha,\beta}\; {{R\left( {\alpha,\beta} \right)}{M^{j}\left( {f_{\alpha}^{(1)},f_{\beta}^{(2)}} \right)}}}},{j = 1},2} & (34) \end{matrix}$

This can be called the value of the two-node competition; a zero sum. The zero sums, for any round where

=(ρ⁽¹⁾,ρ⁽²⁾); is:

P(f ⁽¹⁾ ,f ⁽²⁾;1,1)=P(f ⁽¹⁾ ,f ⁽²⁾;0,0)=0  (35)

But if by p_(ij) I define the probability that Σ¹ loses when the range (ρ_(i) ⁽¹⁾, ρ_(j) ⁽²⁾) is employed and by q_(ij)=1−p_(ij) the probability that Σ¹ wins for the same range, then, for this set of ranges, the expectation of winning g_(ij) for Σ¹ is:

g _(ij) q _(ij) −p _(ij)  (36)

And hence, the rectangular matrix [g_(ij)], i=1, . . . , ρ₁; and j=1, . . . , ρ₂; is the competition matrix.

Consider the competition between Σ¹ and another cluster of nodes that has a fixed range mix (note that the second cluster resembles the set of stationary roadside units). If the second cluster employs another learning automaton, it is fixed range mix would vary slightly over time as its reinforcement scheme alters the state probability. Also note that the alteration would reverse, and then keep oscillating, the reason is that the technique enforces discrete disconnected steps on φ to maintain the finite state of the system.) Let it use the jth range with probability:

$\begin{matrix} {m_{j},{{\sum\limits_{j = 1}^{r_{2}}\; m_{j}} = 1}} & (37) \end{matrix}$

The m_(j) do not depend upon the range chosen by Σ¹. If at the instance t the automaton Σ¹ chooses ρ_(i) ⁽¹⁾, then the probabilities q_(i) of a win by Σ¹ and p_(i) of a loss by Σ¹ are given by:

$\begin{matrix} {{q_{i} = {\sum\limits_{j = 1}^{r_{2}}\; q_{ij}}},m_{j},{p_{i} = {\sum\limits_{j = 1}^{r_{2}}\; p_{ij}}},m_{j}} & (38) \end{matrix}$

Hence, the expected winning for the ith range can be written as:

g _(i) =q _(i) −p _(i)=1−2p _(i)  (39)

Which is identical to Σ¹ being placed in the random medium:

C=C(p ₁ , . . . ,p _(r1)).

Since Σ¹ is a finite deterministic automaton, linear resolutions may apply and can be easily driven. The following result can be shown to be valid. This automaton has ρ₁ valid ranges corresponding to ρ₁ possible branching of the state transition diagram and n states in each branch. If one or more g_(i)≧0, W would represent the asymptotic expected winning of Σ¹ when n→∞ as illustrated here:

W=max(g ₁ , . . . ,g _(r) ₁ )  (40)

If all other nodes obtained optimal ranging estimate, then the value W would be

$\begin{matrix} {W = {\max\limits_{i}\mspace{11mu} {\min\limits_{j}{\sum\limits_{j}\; {g_{ij}m_{j}}}}}} & (41) \end{matrix}$

In such a case, the competition process is deterministic and falls into the same situation of fixed nodes. The cluster of nodes would keep oscillating due to the technique that enforces discrete disconnected steps, which are imposed to maintain the finite state of the system. But that wouldn't hurt the optimality as the competition would oscillate around the optimum, with the small margin selected to convert φ to a discrete set. The smaller steps chosen for φ, the closer (but slower) to optimality. Similarly, larger steps chosen for φ, the wider (but faster) from optimality. One can envision a mixed approach with varying steps, but communicating and synchronizing the change in φ steps would burden the system in practical sense, and would slow the optimization anyway.

Alternatively, if g_(i) is negative, then:

$\begin{matrix} {W = \frac{r_{1}}{\sum\frac{1}{g_{i}}}} & (42) \end{matrix}$

In other words, W would equal to the harmonic mean of g_(i). And the performance would not be optimal following the sense of the competition theory. But equation 41 is valid consistently only when just one g_(i)≧0. When more than one of g_(i)≧0 exists; the probability of suboptimal performance increases consistently. Hence, equation 41 will never be reached, which counters the non-optimal argument.

Comparing Different Automaton

Let us now compare this deterministic automaton's performance with that of a stochastic automaton described earlier. Only the square law reinforcement following the reinforcement scheme described is considered. Then, I let Σ¹ have ρ₁ states and ρ₁ ranges. In state φ_(i), Σ¹ uses the pure range ρ_(i). On the occurrence of this range, the competition environment (i.e. other nodes) uses the jth range with probability m_(j). By following the same arguments given previously, one obtains from the competition matrix G the unit loss probability p_(i) and the unit win probability q_(i). The square law automation performs optimally if any p_(o)<½, in other words, if Π*_(i)=1 and Π*_(j)=0, where j≠i, and the asterisk characterizing the final values. The condition p_(i)<½ means that g_(i)>0. Under such a condition:

W=max(g ₁ , . . . ,g _(r) ₁ )  (43)

Just as described before in equation 40. Therefore, this stochastic finite state automaton has a performance equivalent to the deterministic automaton only when the later has an infinitely large number of states. Again, if the m_(j) are such that they constitute an optimal ranging (using the same multi-lateration approach), W becomes the winning value of the competition as in equation 41.

When the condition g_(i)>0 is not satisfied for both the deterministic and stochastic automata, they play for the harmonic mean of the g_(i). This follows from equation 41 for case of stochastic automaton. In this case, W ends up between the upper and lower values of the competition,

$\begin{matrix} {{\max\limits_{i}\mspace{11mu} {\min\limits_{j}\mspace{11mu} g_{ij}}} \leq W \leq {\min\limits_{j}\mspace{11mu} {\max\limits_{i}\mspace{11mu} g_{ij}}}} & (44) \end{matrix}$

Next two automata zero sum competitions can be considered, that can prove that for deterministic automata of, somewhat, different construction from the linear one and infinitely large number of states. The statement here is modified to exclude the possibility of sub-optimality which may or may not be taken into account. Further proofs of sub-optimality in this type of deterministic automata have been demonstrated by computer experiments and simulations. Sub-optimality can be measured by checking if the following conditions apply:

-   -   1) If the matrix G has one and only one row of positive (column         of negative) elements, then W is the harmonic mean of the         elements of this row (column).     -   2) If G is such that there is no row of all positive elements or         no column of all negative elements, then, W=0.

In exploring competitions between two stochastic automata, Σ¹ and Σ², each with square law reinforcement. Let the state probabilities of Σ¹ be denoted by Π_(i) ⁽¹⁾, i=1, . . . , ρ₁, and Π_(j) ⁽²⁾, j=1, . . . , ρ₂. Final values will be indicated by asterisk. If the probability of unit loss for Σ^(α) is denoted by p_(i) ^((α)), with a range ρ_(i) ^((α)), then:

$\begin{matrix} {p_{i}^{(1)} = {\sum\limits_{j}\; {p_{ij}\Pi_{j}^{(2)}}}} & (45) \end{matrix}$

Where p_(ij)=Pr (unit win for Σ¹ with F={ρ_(i) ⁽¹⁾,ρ_(j) ⁽²⁾)}). And similarly:

$\begin{matrix} {p_{j}^{(2)} = {\sum\limits_{i}\; {\left( {1 - p_{ij}} \right)\Pi_{i}^{(1)}}}} & (46) \end{matrix}$

Like before, the competition matrix G=[g_(ij)] can be used to obtain:

$\begin{matrix} {{p_{ij} = \frac{1 - g_{ij}}{2}},{q_{ij} = \frac{1 + g_{ij}}{2}}} & (47) \end{matrix}$

From equations 45 to 47

$\begin{matrix} {{p_{i}^{(1)} = {\sum\limits_{j}\; \frac{1 - g_{ij}}{2}}},{\Pi_{j}^{(2)} = {\frac{1}{2}\left( {1 - {\sum\limits_{j}\; {g_{ij}\Pi_{j}^{(2)}}}} \right)}}} & (48) \end{matrix}$

Assume that there is a final stochastically stable condition and that the automata have settled in that condition. Now, to determine under what condition p_(i) ^((1)*)<½, which is the condition for optimality in the square law nonlinear reinforcement, from equation48, an equivalent condition is obtained:

$\begin{matrix} {{\sum\limits_{j}\; {g_{ij}\Pi_{j}^{{(2)}^{*}}}} > 0} & (49) \end{matrix}$

Furthermore, it is required that Σ¹ to be optimal for all values of Π(2)*, i.e. whatever be the ranges chosen by Σ². Then equation 49 becomes:

g _(ij)>0, j=1, 2, . . . , r ₂  (50)

Once equation 50 is satisfied, it automatically results in Π^((1)*) no having the ith component unity and the others zero. That's to say Σ¹ is optimal irrespective of the selected ranges in Σ². Equation 50 is equivalent to the existence of an all positive row in the matrix [g_(ij)]. Since [g_(ij)] corresponds to the winnings of Σ¹, modification of equation 50 in order to make Σ² optimal results in the condition:

g _(ij)<0, i=1,2, . . . ,r ₁  (51)

This is equivalent to a column with all negative elements. Further, suppose that equation 51 is satisfied for some value of i. Then, equation 50 cannot be satisfied, i.e., both Σ¹ and Σ² cannot be simultaneously asymptotically optimal. Now if Σ¹ is optimal, then Σ² has a final state probability distribution given to a great degree of accuracy by:

$\begin{matrix} {\Pi_{j}^{{(2)}^{*}} = \frac{\frac{1}{{2p_{j}^{(2)}} - 1}}{\left( {\sum\limits_{j}\; \frac{1}{{2p_{j}^{(2)}} - 1}} \right)}} & (52) \end{matrix}$

Then, by reapplying the optimality condition of Σ¹ from equation 46:

$\begin{matrix} {p_{j}^{{(2)}^{*}} = {{1 - p_{ij}} = \frac{1 + g_{ij}}{2}}} & (53) \end{matrix}$

Therefore, equation 52 can be rewritten as:

$\begin{matrix} {\Pi_{j}^{{(2)}^{*}} = \frac{\frac{1}{g_{ij}}}{\sum\limits_{j}\; \frac{1}{g_{ij}}}} & (54) \end{matrix}$

Equation 54 can be used to drive the average winnings of Σ¹ as follows:

$\begin{matrix} \begin{matrix} {W = \frac{r_{1}}{\sum\limits_{j}\; \frac{1}{g_{ij}}}} \\ {= {{the}\mspace{14mu} {harmonic}\mspace{14mu} {mean}\mspace{14mu} {of}\mspace{14mu} {the}\mspace{14mu} {jth}\mspace{14mu} {{row}.}}} \end{matrix} & (55) \end{matrix}$

Now interchanging the two automata Σ¹ and Σ², one can easily establish that:

if g _(ij)<0, i=1,2, . . . r ₁  (56)

Then:

$\begin{matrix} \begin{matrix} {W = \frac{r_{2}}{\sum\limits_{i}\; \frac{1}{g_{ij}}}} \\ {= {{the}\mspace{14mu} {harmonic}\mspace{14mu} {mean}\mspace{14mu} {of}\mspace{14mu} {the}\mspace{14mu} {jth}\mspace{14mu} {{row}.}}} \end{matrix} & (57) \end{matrix}$

A view on how the entire process is performed is elaborated in FIG. 6 where the received position announcements are used along with the GNSS estimated position to calculate the accurate position and broadcast it such as in a dedicated short-range communications (DSRC) vehicular environment. The system 600 may be operable in a vehicle where DSRC communication or other inter-vehicle communication system. The vehicle receives position data via a GNSS system using for example GPS from a GNSS network 630 through a GNSS receiver 610. In addition, the DSRC receiver 612 provide position announcements from vehicles 640 that are in proximity of the vehicle containing system 600. As discussed above the received position announcements are processed to refine the accuracy of the vehicles current position. This can be performed by processing the position announcements by one or more processors 622 containing instructions in memory 620 for processing the data. Additional ranging data may be provided by ranging systems such as Lidar (Light detection and ranging), radar, or any wireless or optical ranging technology (not shown) which would enable correlation between the received position announcements and a vehicle position. The memory may comprise modules such as time synchronization module 614 for correlating the position announcements and GPS position data, a fuzzy clustering module 616 for performing clustering functions on the position data, and a stochastic automata model module 618 for determining a more accurate position based upon the received position announcements. The improved position determination can then be transmitted by a DSRC transmitter 624 to other vehicles.

FIG. 7 shows system for performing cooperative stochastic positioning in a mobile environment, for example a closed environment such as a store or warehouse where GNSS systems may not work or provide the required resolution. The system 700 may be operable in any environment that absolute or relative position information is required such as in inventory or warehouse situations using radio frequency identification (RFID) technologies, mobile device positioning or any relative object position applications within a defined coordinate system. The object receives position data via position data receiver 710. The position data may include only relative object position data or may include reference position data based upon fixed position points in the immediate area. In addition, GPS position data may also be utilized if available. As discussed above the received position announcements are processed to refine the accuracy of the objects current position. This can be performed by processing the position announcements by one or more processors 720 containing instructions in memory 718 for processing the data. The memory may comprise modules such as time synchronization module 712 for correlating the position announcements, a fuzzy clustering module 714 for performing clustering functions on the position data, and a stochastic automata model module 716 for determining a more accurate position based upon the received position announcements. The improved position determination can then be transmitted by a position data transmitter 722 to other objects. Position announcements may be transmitted or received by wireless technologies such as Bluetooth, Zigbee, IEEE 802.11, 802.16, cellular or mobile wireless technologies.

Specific embodiments have been shown and described herein. However, modifications and variations may occur to those skilled in the art. All such modifications and variations are believed to be within the scope and sphere of the present disclosure. 

1. A method for cooperative stochastic positioning in a mobile environment executed by a processor in an object of interest, the method comprising: receiving position announcements transmitted wirelessly from a plurality of objects, the position announcements providing the current position data of the respective object in relation to a common coordinate system; discretizing the received position announcements to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data; performing clustering of the data groupings to determine which cluster datasets from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the object of interest; applying stochastic automata model to selected cluster dataset to evaluate relative cluster weights in order to determine the relative position of the object of interest; and updating accuracy of a current position of the object of interest based upon determined relative position.
 2. The method of claim 1 wherein discretizing the received position data further comprises applying a channel propagation model on the received object position data to determine range data associated with the respective received object position data.
 3. The method of claim 1 wherein discretizing the received position data further comprises receiving range data associated with the received object position data associated with the respective received object position data.
 4. The method of claim 1 wherein the performing clustering further comprises: applying classifiers to the determined clusters based upon selected criteria; determining if the selected cluster dataset provides object position data accuracy within a defined position tolerance; and dividing the cluster dataset if the determined accuracy is less than the defined tolerances.
 5. The method of claim 1 wherein the performing clustering further comprises: selecting clustering criteria to be applied to the received positioning announcements; selecting a cluster dataset from the data groupings; minimizing error of the selected cluster dataset using selected criteria; determining average normalized intergroup scatter between position data provided in the cluster dataset; determining a normalized scatter matrix of the position data provided in the cluster dataset; determining if criteria should be adjusted based upon normalized scatter matrix; determining if cluster dataset should be combined with previously analysed cluster based upon a defined criteria accuracy; determining if the cluster dataset provides object position data accuracy within defined position tolerances; and dividing the cluster dataset if the determined accuracy is less than the defined tolerances to perform clustering on each of the divided clusters.
 6. The method of claim 1 wherein applying stochastic automata model to selected clusters to determine the relative position further comprises: calculating a weighted multi-lateration position using the provided cluster datasets; applying weighted ranges and a set of multi-lateration to determine an accurate position from the position data in the cluster dataset; validate accuracy of obtained position relative to current object of interest position wherein least weighted data is removed from the dataset when dataset positions accuracy is less than required accuracy threshold and current position data is updated with obtained position.
 7. The method of claim 1 wherein the clustering criteria is a function of position announcement data properties either determinable by the received announcement data or associated with previously determined characteristics associated with the object providing the position announcement.
 8. The method of claim 7 wherein the clustering criteria is selected from the group comprising: least-mean-square (LMS), arrival delay and position announcement data confidence.
 9. The method of claim 1 further comprising receiving positioning information at the object of interest from a global positioning system (GPS) and wherein receiving the positioning announcements transmitted from the plurality of objects is via a short range inter-object communication system and the relative object positioning is determined relative to a fixed coordinate system. 10-11. (canceled)
 12. The method of claim 9 wherein the wireless technology is selected from the group comprising: Bluetooth, Zigbee, IEEE 802.11, 802.16, and cellular based wireless communication technologies.
 13. A system for cooperative stochastic positioning in a mobile environment the system comprising: a first receiver for receiving positioning announcements from a plurality of objects, the position announcements providing the current position data of the respective object in relation to a common coordinate system; a processor coupled to the first receiver; and a memory comprising instructions for execution by the processor, the instructions comprising: discretizing the received position announcements to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data; performing clustering of the data groupings to determine which cluster datasets from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the object of interest; applying stochastic automata model to selected cluster dataset to evaluate relative cluster weights in order to determine the relative position of an object of interest; and updating accuracy of a current position of the object of interest based upon determined relative position.
 14. The system of claim 13 further comprising a second receiver for receiving positioning data from a global navigation satellite system for determining a position of the object of interest relative to a defined global coordinate system for the current position.
 15. The system of claim 14 wherein the first receiver is a short range wireless communication system for communication between objects in an immediate are near the object of interest.
 16. The system of claim 15 wherein the short range wireless communication system is a dedicated short-range communications (DSRC) for vehicular applications.
 17. The system of claim 13 further wherein receiving the positioning announcements transmitted from the plurality of objects is a dedicated short range wireless technology and the relative object positioning is determined relative to a fixed coordinate system.
 18. The system of claim 17 wherein the wireless technology is selected from the group comprising: Bluetooth, Zigbee, IEEE 802.11, 802.16, and cellular based wireless communication technologies.
 19. The system of claim 13 wherein discretizing the received position data further comprises applying a channel propagation model on the received object position data to determine range data associated with the respective received object position data.
 20. The system of claim 13 wherein discretizing the received position data further comprises receiving range data associated with the received object position data associated with the respective received object position data.
 21. The system of claim 13 wherein the performing clustering further comprises: applying classifiers to the determined clusters based upon selected criteria; determining if the selected cluster dataset provides object position data accuracy within a defined position tolerance; and dividing the cluster dataset if the determined accuracy is less than the defined tolerances.
 22. The system of claim 13 wherein the performing clustering further comprises: selecting clustering criteria to be applied to the received positioning announcements; selecting a cluster dataset from the data groupings; minimizing error of the selected cluster dataset using selected criteria; determining average normalized intergroup scatter between position data provided in the cluster dataset; determining a normalized scatter matrix of the position data provided in the cluster dataset; determining if criteria should be adjusted based upon normalized scatter matrix; determining if cluster dataset should be combined with previously analysed cluster based upon a defined criteria accuracy; determining if the cluster dataset provides object position data accuracy within defined position tolerances; and dividing the cluster dataset if the determined accuracy is less than the defined tolerances to perform clustering on each of the divided clusters.
 23. The system of claim 13 wherein applying stochastic automata model to selected clusters to determine the relative position further comprises: calculating a weighted multi-lateration position using the provided cluster datasets; applying weighted ranges and a set of multi-lateration to determine an accurate position from the position data in the cluster dataset; validate accuracy of obtained position relative to current object of interest position wherein least weighted data is removed from the dataset when dataset positions accuracy is less than required accuracy threshold and current position data is updated with obtained position.
 24. The system of claim 13 wherein the clustering criteria is a function of position announcement data properties either determinable by the received announcement data or associated with previously determined characteristics associated with the object providing the position announcement.
 25. The system of claim 24 wherein the clustering criteria is selected from the group comprising: least-mean-square (LMS), arrival delay and position announcement data confidence.
 26. A computer readable memory providing instructions for performing cooperative stochastic positioning in a mobile environment, the when executed by a processor in an object of interest performing the method comprising: receiving wirelessly, position announcements transmitted from a plurality of objects, the position announcements providing the current position data of the respective object in relation to a common coordinate system; discretizing the received position announcements to obtain a plurality of data groupings based upon a clustering criteria applied to the received position data; performing clustering of the data groupings to determine which cluster datasets from the data groupings provide sufficient and consistent position accuracy to determine a relative position of the object of interest; applying stochastic automata model to selected cluster dataset position to evaluate relative cluster weights in order to determine the relative position of the object of interest; and updating accuracy of a current position of the object of interest based upon determined relative position. 